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. Abstract. I present a selection of recent lattice results in flavourdynamics, including the status of 

the calculation of quark masses and a variety of weak matrix elements relevant for the determination 
of CKM matrix elements. Recent improvements in the momentum resolution of lattice computations 
and progress towards precise computations of K — ► %n decay amplitudes are also reviewed. 

o : 

INTRODUCTION 

7— I ' 

■^j- \ One of the main approaches to testing the Standard Model of Particle Physics and search- 

ing for signatures of new physics is to study a large number of physical processes to 
obtain information about the unitarity triangle and to check its consistency. The preci- 

O ! si° n with which this check can be accomplished is limited by non-perturbative QCD 

effects and lattice QCD provides the opportunity to quantify these effects without model 
assumptions. Of course, lattice computations themselves have a number of sources of 
systematic uncertainty, and much of our current effort is being devoted to reducing and 
controlling these errors. In this talk I briefly discuss the evaluation of quark masses and 
weak matrix elements using lattice simulations. 
r-| ! For most lattice calculations of physical quantities, the principal source of systematic 

uncertainty is the chiral extrapolation, i.e. the extrapolation of results obtained with 
■ unphysically large u and d quark masses. Ideally we would like to perform computations 

with 140 MeV pions and hence with m q /m s of about 1/25 (where m q (m s ) is the average 
light quark mass (strange quark mass)). In practice values m q /m s > 111 are fairly 
typical, so that the MILC Collaboration's simulation with m q /m s ~ 1/8 is particularly 
impressive IH and provides a challenge to the rest of the community to reach similarly 
low masses. Its configurations have been widely used to determine physical quantities 
with small quoted errors. 

The MILC collaboration uses the staggered formulation of lattice fermions and for a 
variety of reasons it is very important to verify the results using other formulations. With 
staggered fermions each meson comes in 16 tastes and the unphysical ones are removed 
by taking the fourth root of the fermion determinant. Although there is no demonstration 
that this procedure is wrong, there is also no proof that it correctly yields QCD in the 
continuum limit 01 . The presence of unphysical tastes leads to many parameters to be fit 
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in staggered chiral perturbation theory (typically many tens of parameters) and to date 
the renormalization has only been performed using perturbation theory. It is therefore 
pleasing to observe that the challenge of reaching lower masses is being taken up by 
groups using other formulations of lattice fermions (see e.g. ref. . 

In this talk I will discuss a selection of issues and results in lattice flavourdynamics. 
I start by describing some new thoughts on improving the momentum resolution in 
simulations, by varying the boundary conditions on the quark fields. I then review the 
status of lattice calculations of quark masses, decays (for which computations have 
only recently began) and Bk- This is followed by a discussion of some of the key issues 
in the computation of K — > %% decays and in heavy-quark physics. 



Improving the Momentum Resolution on the Lattice 

Numerical simulations of lattice QCD are necessarily performed on a finite spatial 
volume, V = I? . Providing that V is sufficiently large, we are free to choose any 
consistent boundary conditions for the fields (j)(x,t), and it is conventional to use periodic 
boundary conditions, <j>(xj + L) = (jc/) (i = 1,2 or 3). This implies that components 
of momenta are quantized to take integer values of 2%/L. Taking a typical example 
of a lattice with 24 points in each spatial direction, L = 24a, with a lattice spacing 
a = 0.1 fm so that aT x ~ 2GeV, we have 2n/L = .52GeV. The available momenta 
for phenomenological studies (e.g. in the evaluation of form-factors) are therefore very 
limited, with the allowed values of each component pj separated by about 1/2 GeV. The 
momentum resolution in such simulations is very poor. 

Bedaque [4] has advocated the use of twisted boundary conditions for the quark fields 
q(x) e.g. 

■a 2n 0j 

q(xi+L) = e l q(xi) with momentum spectrum Pi = ni 1 , (1) 

L L 

with integer n,;. Modifying the boundary conditions changes the finite-volume effects, 
however, for quantities which do not involve Final State Interactions (e.g. hadronic 
masses, decay constants, form-factors) these errors remain exponentially small also with 
twisted boundary conditions [5]. Since we usually neglect such errors when using peri- 
odic boundary conditions, we can use twisted boundary conditions with the same preci- 
sion. Moreover the finite- volume errors are also exponentially small for partially twisted 
boundary conditions in which the sea quarks satisfy periodic boundary conditions but 
the valence quarks satisfy twisted boundary conditions o,|g]. This is of significant prac- 
tical importance, implying that we do not need to generate new gluon configurations for 
every choice of twisting angle { 0,- } . 

The use of partially twisted boundary conditions opens up many interesting phe- 
nomenological applications, solving the problem of poor momentum resolution. It also 
appears to work numerically. Consider for example, the plots in fig.E obtained using an 
unquenched (2 flavours of sea quarks) UKQCD simulation on a 16 3 x 32 lattice, with a 
spacing of about 0. 1 fm. The plots correspond to a value for the light-quark masses for 
which m n /m p = 0.7 0]. The lower (upper) left-hand plot shows the energy of the % 
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FIGURE 1. Plots of the Dispersion Relation (left) and Decay Constants (right) as a function of the 
momentum p of the mesons. In both cases the top (bottom) plot corresponds to the p-meson (7T-meson). 



(p) as a function of the momentum of the meson, and the right-hand plot shows the bare 
values of the leptonic decay constants f n and f p . The x-axis denotes (\p\L) 2 . The re- 
sults are beautifully consistent with expectations (particularly for pL < 2% where lattice 
artifacts are small); the predicted dispersion relation is satisfied and the extracted decay 
constants are independent of the momenta. Using periodic boundary conditions only the 
results at values of p indicated by the dashed lines are accessible. With partially twisted 
boundary conditions all momenta are reachable. 



QUARK MASSES 

Quark Masses are fundamental parameters of the Standard Model, but unlike leptons, 
quarks are confined inside hadrons and are not observed as physical particles. Quark 
masses therefore cannot be measured directly, but have to be obtained indirectly through 
their influence on hadronic quantities and this frequently involves non-perturbative QCD 
effects. Lattice simulations prove to be very useful in the determination of quark masses; 
particularly for the light quarks (u, d and s) for which perturbation theory is inapplicable. 

In order to determine the quark masses we compute a convenient and appropriate set 
of physical quantities (frequently a set of hadronic masses) and vary the input masses 
until the computed values correctly reproduce the set of physical quantities being used 
for calibration. In this way we obtain the physical values of the bare quark masses, from 
which by using perturbation theory, or preferably non-perturbative renormalization, the 
results in standard continuum renormalization schemes can be determined. 

My current best estimates for the values of the quark masses as determined from 
lattice simulations are presented in table [T] 

The relatively large error on the mass of the charm quark is a reflection of the fact 
that the most detailed study to date was performed in the quenched approximation [|8|], 
whose authors find m c = 1.301(34) GeV. I have added a conservative 15% error as an 
estimate of quenching effects. Current and future calculations will be dominated by 
unquenched simulations so that the error will decrease very significantly. Indeed a very 



TABLE 1. My summary of the status of lattice determinations of quark masses (in the MS renormal- 
ization scheme). For m = (m u + m c /)/2 and m s the results are presented at 2GeV and for m c , nib, the 
results are presented at the mass itself (in = m(m)) . For comparison the values quoted by the PDG in 
2004, using or excluding lattice simulations, are also presented. 



Flavour 


Best Lattice Values 


PDG 2004 (Lattice) 


PDG 2004 (Non-Lattice) 


»i(2GeV) 

m s (2GeV) 
m c 
m b 


(3.8±0.8)MeV 

(95±20)MeV 
(1.26±0.13±0.20)GeV 
(4.2±0.1±0.1)GeV 


(4.2±1.0)MeV 

(105±25)MeV 
(1. 30 ± 0.03 ± 0.20) GeV 
(4.26±0.15±0.15)GeV 


( 1 .5 < m u (2 GeV) < 5) Me V 
(5 <m d (2GeV) <9)MeV 
80 - 155 Me V 
l-1.4GeV 
4 -4.5 GeV 



recent unquenched calculation finds m c = 1.22(9) GeV 

The relative error on is small because what is actually calculated is mg — mi,. 
The calculations are performed in the Heavy Quark Effective Theory and the major 
source of systematic error is the subtraction of 0(1/ (oAqq))) terms. Using stochastic 
perturbation theory, Di Renzo and Scorzato have performed this calculation to 3-loop 
order ifioll . The second error on fht, in table[TJis my conservative estimate of the fact that 
the simulations have been performed with two flavours of sea quarks. 



SELECTED TOPICS IN KAON PHYSICS 

Decays 

A new area of investigation for lattice simulations is the evaluation of non-perturbative 
QCD effects in K — > n£vg decays, from which the CKM matrix element V us can be 
determined. The QCD contribution to the amplitude is contained in two invariant form- 
factors f°(q 2 ) and f + (q 2 ) defined by 

( Tt(p K ) \s Yll u | K(p K ) ) = f°(q 2 ) K 2 % +f+(q 2 ) 

q L 

where q = Pk~ Pk- (Parity Invariance implies that only the vector current from the V — A 
charged current contributes to the decay.) A useful reference value for /+ (0) comes from 
the 20-year old prediction of Leutwyler and Roos, f+(0) = 1 + f 2 + f 4 + ■ ■ ■ = 0.961(8) 
where f n = 0{M\ 7trj ).f2 = —0.023 is well determined, whereas the higher order terms 
in the chiral expansion require model assumptions. 

To be useful in extracting V us from experimental measurements we need to be able 
to evaluate /°(0) = / + (0) to better than about 1% precision. This would seem to be 
impossible until one notes that it is possible to compute 1 — / + (0), so that an error of 
1% on / + (0) is actually an error of 0(25%) on 1 — / + (0). The calculation follows a 
similar strategy to that proposed in ref. Ill ill for the form-factors ofB^D semileptonic 
decays (which in the heavy quark limit are also close to 1), starting with a computation 
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FIGURE 2. A compilation of recent quenched results for B^ s (2 Ge V) . 



of double ratios such as 

(7t\sYol\K)(K\lY Q s\7t) _ r Q 2 -.2 (m K + m K ) 2 
{n\ly Q l\K)(K\sy Q s\K) [J ^ maxJJ 4m K m K ' ^ ) 

where all the mesons are at rest and g max = (Mk—M^) 2 . 

Following a quenched calculation by the SPQR collaboration last year lfl3ll . in which 
the strategy for determining the form-factors was presented, there have been 3 very 
recent unquenched (albeit largely preliminary) results: 

RBC[HJ / + (0) = 0.955 (12) 
JLQCD [15J /+(0) = 0.952 (6) 
FNAL/MILC/HPQCD □£] /+ (0) = 0.962 (6) (9) 

in good agreement with the result of Leutwyler and Roos lfl2ll . 



B K 

Bk, the parameter which contains the non-perturbative QCD effects in K° — K° mix- 
ing, has been computed in lattice simulations by many groups. It is defined by 

(K° | (sy^(l — y 5 )^) -fWlK ) = \Mlf K B K . (3) 

Bk depends on the renormalization scheme and scale and is conventionally given in the 
NDR, MS scheme at \x = 2 GeV or as the RGI parameter B% (the relation between the 

two is B K ~ 1.4 5M S (2GeV)). 

A compilation of recent results for Bk obtained in the quenched approximation is 
presented in fig.El From such results recent reviewers have summarised the status of 



quenched calculations as: 



B% s (2GeV) =0.58(4) O and flf 5 (2GeV) = 0.58(3) GjJ . (4) 

The dashed lines in fig. correspond to Bk = .58(3), which I am happy to take as the 
current best estimate. 

The challenge now is to obtain reliable unquenched results; such computations are un- 
derway by several groups but so far the results are very preliminary. We will have to wait 
a year or two for precise results, but I mention in passing C.Dawson's guesstimate 111 811 
(stressing that it is only a guesstimate), based on a comparison of quenched and un- 
quenched results at similar masses and lattice spacings, of 5^ 5 (2GeV) =0.58(3)(6). 

K — > kk Decays 

A quantitative understanding of non-perturbative effects in K — > kk decays will be an 
important future milestone for lattice QCD. Two particularly interesting challenges are: 

i) an understanding of the empirical AI = 1/2 rule, which states that the amplitude for 
decays in which the two-pion final state has isospin 1=0 is larger by a factor of about 22 
than that in which the final state has 1 = 2; 

ii) a calculation of e'/e, whose experimental measurement with a non-zero value, 
(17.2 ±1.8) x 10~ 4 , was the first observation of direct CP- violation. 

The two challenges require the computation of the matrix elements of the AS = 1 oper- 
ators which appear in the effective Weak Hamiltonian. 

About 4 years ago, two collaborations published some very interesting quenched 
results for these quantities: 



Collaboration(s) 


Re A /Re A 2 


e'/e 


RBC [12J 
CP-PACS[20J 
Experiments 


25.3 ±1.8 
9-12 
22.2 


-(4.0±2.3) x 10" 4 

(-7--2)xl0~ 4 
(17.2±1.8) x 10" 4 



Both collaborations obtain a considerable octet enhancement (significantly driven how- 
ever, by the chiral extrapolation) and e' J £ with the wrong sign. A particularly impressive 
feature of these calculations was that the collaborations were able to perform the sub- 
traction of the unphysical terms which diverge as powers of the ultra-violet cut-off (a , 
where a is the lattice spacing). The results are very interesting and will provide valuable 
benchmarks for future calculations, however the limitations of the calculations should be 
noted, in particular the use of chiral perturbation theory CfPT) only at lowest order. This 
has the practical advantage that K — > kk matrix elements do not have to be evaluated 
directly, it is sufficient at lowest order to study the mass dependence of the matrix ele- 
ments (M\ &i \M) and (0| &i\M), where M is a pseudoscalar meson and the G\ are the 
AS = 1 operators appearing in the effective Hamiltonian, to determine the low-energy 
constants and hence the amplitudes. It is not very easy to estimate the errors due to this 
approximation, but they should be at least of 0(mj, / Aq^) . Since for e'/e the dom- 
inant contributions appear to be from the QCD and electroweak penguin operators &(, 



and 6%, which are comparable in magnitude but come with opposite signs, it is not to- 
tally surprising that the prediction for e'/ e at lowest order in #PT has the wrong sign. 
It should also be noted that in the simulations described in ref. H19L the light quarks 
masses were large (the pions were heavier than about 400 MeV) and so one can question 
the validity of ^PT in the range of masses used (about 400-800 MeV). 

To improve the precision, apart from performing unquenched simulations and reduc- 
ing the masses of the light quarks, one needs to go beyond lowest order ^PT (for exam- 
ple by going to NLO l!2lL 12211 ) and, in general, this requires the evaluation of K — > %n 
matrix elements and not just M — > M ones. The treatment of two-hadron states in lattice 
computations has a new set of theoretical issues, most notably the fact that the finite- 
volume effects decrease only as powers of the volume and not exponentially. Starting 
with the pioneering work of Liischer [23], the theory of finite- volume effects for two- 
hadron states in the elastic regime is now fully understood, both in the centre-of-mass 
and moving frames, 112311 - B28J and I will now briefly discuss this. 

Consider the two-hadron correlation function represented by the diagram 




where the shaded circles represent two-particle irreducible contributions in the s- 
channel. For simplicity let us take the two-hadron system to be in the centre-of mass 
frame and assume that only the s-wave phase-shift is significant (the discussion can be 
extended to include higher partial waves). Consider the loop integration/summation over 
p (see the figure). Performing the po integration by contours, we obtain a summation 
over the spatial momenta of the form: 

p 1 

where the relative momentum k is related to the energy by E 2 — 4(m 2 + k 2 ) , the function 
f(p 2 ) is non-singular and (for periodic boundary conditions) the summation is over 
momenta p = (2n/L)n where n is a vector of integers. In infinite volume the summation 
in eq. © is replaced by an integral and it is the difference between the summation 
and integration which gives the finite- volume corrections. The relation between finite- 
volume sums and infinite-volume integrals is the Poisson Summation Formula, which 
(in 1-dimension) is: 

p l=— °° 

If the function g(p) is non-singular, the oscillating factors on the right-hand side ensures 
that only the term with 1 = contributes, up to terms which vanish exponentially with 
L. The summand in eq. © on the other hand is singular (there is a pole at p 2 = k 2 ) 
and this is the reason why the finite- volume corrections only decrease as powers of L. 
The detailed derivation of the formulae for the finite- volume corrections can be found 
in refs. ll23ll - ll28ll and is beyond the scope of this talk. The results hold not only for 
K — ► %% decays, but also for % - nucleon and nucleon-nucleon systems. 



For decays in which the two-pions have isospin 2, we now have all the necessary 
techniques to calculate the matrix elements with good precision and such computations 
are underway. For decays into two-pion states with isospin there are also no barriers 
in principle. However, in this case, purely gluonic intermediate states contribute and we 
need to learn how to calculate the corresponding disconnected diagrams with sufficient 
precision. In addition the subtraction of power- like ultraviolet divergences requires large 
datasets (as demonstrated in refs. [ 19, 20] in quenched QCD). For these reasons it will 
take a longer time for some of the AI = 1/2 matrix elements to be computed than 
AI = 3/2 ones. 

HEAVY QUARK PHYSICS 

Lattice simulations are playing an important role in the determination of physical quanti- 
ties in heavy quark physics including decay constants (fB,fB s ,fD,fD s ), the 5-parameters 
of B — B mixing (from which the CKM matrix elements V t d and V ts can be determined), 
form-factors of semileptonic decays (which give V c b and V u b), the gBB*n coupling con- 
stant of heavy-meson chiral perturbation theory and the lifetimes of beauty hadrons. 

The typical lattice spacing in current simulations a ~ 0. 1 fm is larger than the Comp- 
ton wavelength of the £>-quark and comparable to that of the c-quark. The simulations 
are therefore generally performed using effective theories, such as the Heavy Quark Ef- 
fective Theory or Non-Relativistic QCD. Another interesting approach was proposed by 
the Fermilab group 12911 . in which the action is improved to the extent that, in principle 
at least, artefacts of 0((niQa)") are eliminated for all n, where Mq is the mass of the 
heavy quark Q. Determining the coefficients of the operators in these actions requires 
matching with QCD, and this matching is almost always performed using perturbation 
theory (most often at one-loop order). This is a significant source of uncertainty and 
provides the motivation for attempts to develop non-perturbative matching techniques. 

I only have time here to consider very briefly a single topic, semileptonic 5-decays. 
For B — > % decays, the pion's momentum has to be small in order to avoid large lattice 
artefacts, so that q 2 = (ps — Pn) 2 is large (q 2 > 15 GeV 2 or so). There continues to be 
a considerable effort in extrapolating these results over the whole q 2 range. Recently, as 
experimental results begin to be presented in q 2 bins, it has become possible to com- 
bine the lattice results at large q 2 with the binned experimental results and theoretical 
constraints to obtain V u b with good precision [30]. 

As an example I present a recent result, obtained using the MILC gauge field config- 
urations with staggered light quarks and the Fermilab action for the £>-quark 1. 3111 

\V ub \ = 3.48(29) (38) (47) x 10~ 3 . (7) 

I mention that other semileptonic decays of heavy mesons are also being studied, in- 
cluding B -> dW decays (a recent result is \V cb \ =3.9(1)(3) x 10~ 2 0]) and£>^ %,K 
decays. 



SUMMARY AND CONCLUSIONS 



Lattice QCD simulations, in partnership with experiments and theory, play a central 
role in the determination of the fundamental parameters of the Standard Model (e.g. 
quark masses, CKM matrix elements) and in searches for signatures of new physics 
and ultimately perhaps will help to unravel its structure. With the advent of unquenched 
simulations, a major source of uncontrolled systematic uncertainty has been eliminated 
and the main aim now is to control the chiral extrapolation and reduce other systematic 
uncertainties. We continue to extend the range of applicability of lattice simulations to 
more processes and physical quantities. In this talk I have only been able to give a small 
selection of recent results and developments; a more complete set can be found on the 
web-site of the 2005 international symposium on lattice field theory ll32ll . 
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